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1 Introduction 

The Korteweg-de Vries (abbreviated KdV henceforth) equation 

du d 3 u du 
dt dx 3 dx 

for a function u = u(x, t) of position x and time t naturally extends to an infinite sequence 
of partial differential equations of solitonic character |] 

£ = |g»H (. = 1,2,...), (2) 

where each Gk = Gk[u] is some differential polynomial in u with respect to X, 1.6. cL 
polynomial in u and its derivatives 

u„:=|^, (n=l,2,...). (3) 



dx r 



The sequence begins with 



Gi = u, 

G 2 = u 2 + 3u 2 , 

G 3 = u 4 + Wuu 2 + 5uj + 10m 3 , ... (4) 

While usually the KdV hierarchy is defined through a recursive rule for the Gk, we derived 
in our papers [Q, [3| an explicit expression for the latter. The main aim of the present paper 
is to improve this expression by decomposing it in homogeneous parts. Thus, the rest of 
impliciteness in 0, [J is overcome. We derive the fully explicit expression by means of a 
" method of matrix elements" , developed in ||, U , which will be described shortly in section 
2. The same method produces an expression for the powers L k of the one- dimensional 
Schrddinger operator 

d 2 

and gives some knowledge on the differential operators A k in the Lax representation 

r)T 

^- = [A k ,L], (k = 1,2,...) (6) 

of the KdV hierarchy 0. More generally, the method gives, as a side-result, an expression 
for the powers L k of an r-th order Schrddinger-like operator 

QV 

+ u. (7) 



dx 



A particularly convenient formula, in terms of the so called Bell differential polynomials, 
is valid for r = 1. 
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While the paper || treated a scalar potential u = u(x, t) only, the derivation in || and 
in the present paper holds for the matrix case too, where u = u(x,t) is a N x N matrix 
valued function. Notice that the matrix KdV hierarchy (|55|) has solitonic character if u 
is a Hermitian (in particular, a real symmetric) matrix: the papers |6|, [7|, ||, ^ [HJ showed 
that the inverse scattering method for solving Cauchy's problem works well then. 

It is well known (cf., e.g., ]§ g (llj, 0) that the heat kernel K — K{ s; x, x'), i.e. the 
fundamental solution of the heat equation 

(!-i)* = 0, (8) 

to the one dimensional Schrodinger operator @ has an asymptotic expansion as t — > 
of the form: 

K = K(s; x,x') ~ (4 7 r S )- 1 / 2 exp (x £ i ajfc (x, (9) 

The coefficients a k = a,k(x,x') satisfy the recursive differential equations 

1 + -D) a k = La fc _i, (A; =1,2,...), a := 1, (10) 
where we abbreviate 

^:=(x-x')^. (11) 



It is further known 113, 14, H that the right-hand sides of the KdV hierarchy are deter- 



mined by the diagonal values of the heat kernel coefficients: 

1 f2k\ 

G k [u](x) = 2(1. J a k(x,x). (12) 

We will evaluate the recursion system (pTOj ) in order to calculate the diagonal values 
ctk(x, x). In fact, we do more and calculate all Taylor coefficients 

d n 

< n \ a k>--= g^a k (x,x')\ x=xl . (13) 

2 The method of matrix elements 

Let us briefly describe a " method of matrix elements" as developed in |5|, f|, |12| . Let us 
fix a point x' on the real line R and consider a sufficiently small interval around x', say 
I = [x' — R, x' + R] , where R is a sufficiently small positive constant. We consider the 
space C U (I) of analytic functions f = f(x, x') of the real variable x G I and their Taylor 
series 

/(*, x') = J2-(x- xTihf(x, x') . (14) 



n=Q n\ dx n ' 
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We rewrite the Taylor formula ([14]) in the form 

oo 

/=|/>=X>Xn|/>, (15) 



n=0 

where 



\n>:=-(x-x') n , (n = 0,1,2,...), (16) 



(17) 



In a sense, the functions \n > form a basis of the space (7^(1) and the < n| are the 
corresponding dua£ linear functionals < n| : C W (I) — > R. Taylor's theorem can formally 
be expressed as a completeness relation 

oo 

J2\n><n\ = l, (18) 

n=0 

where 1 denotes the identical operator. The duality of the systems < m\ and \n > is 
expressed by 

< m \n >= S mn , (19) 

where 5 mn is the Kronecker symbol. 

Let us consider a linear differential operator of order r with respect to x 6 I 

L = ^u j (a;,x')^— 

with analytic coefficients '). Such an operator can be characterized by its 

matrix elements 

min(r ' n) ( m \ 

< m\L\n >:=< m\(L\n >) = ^ y .j<m — n + j\uj> (20) 



Namely, there holds 



L = 2_, \ m >< m\L\n >< n\, (21) 

n,m>0 



which is just an abbreviation of 

oo 

L\f >= ^2 \m >< m\L\n >< n\f > . (22) 

n,m>0 

Note that for differential operators of order r the matrix elements < m\L\n > vanish for 
n > m + r + 1: 

< m|L|n >= for n > m + r + 1. (23) 



/. G. Avramidi and R. Schimming: Explicit expression for KdV hierarchy 



4 



Therefore, the sum over n in (^T[) is always finite! The infinite sum over m converges 
because of the analyticity of the coefficients Uj of the operator L. This is a special property 
of differential operators and it is not necessarily valid for pseudodifferential operators. 

Formally, fl2~T|) is obtained by application of the identity operator (|T8| ) on both sides. 
By inserting the completeness relation one can show that the matrix elements of 
the product of two operators L X L 2 result from the matrix multiplication of the matrix 
elements of L x and L 2 : 

oo 

< m\L 2 Li\n >= ^ < m\L 2 \l >< l\L x \n > . (24) 
z=o 

Again, from the property ([23|), we see that this sum is finite and there is no problem 
of convergence. The same argument applies to the product L k ■ ■ ■ L 2 Li of operators 
L\, L 2l • • • , L k : 

< m\L k - ■ ■ L 2 L 1 \n> = ^ < m|L fc |n fc _i >< n k „i\L k _i\n k _ 2 > ■ ■ ■ 

n 1 ,n 2 ,...,n k _ 1 >0 

x • ■ ■ < n 2 \L 2 \n\ >< ni\L\\n > . (25) 

Let us note that the method of matrix elements applies also to the more general case 
of smooth functions, i.e. to the space C°°(l). In this case the Taylor series does not, in 
general, converge but is an asymptotic expansion. Thus, all the equations of this section 
hold also for operators with smooth , not necessarily analytic, coefficients, if one replaces 
the sign = by the sign for asymptotic equality ~ in the Taylor expansion. The equations 
between Taylor coefficients remain exactly valid even in this case. 

Let us consider now some examples of linear operators. 

1. We already know that the identical operator 1 has the matrix elements 

< m|l|n >=< m\n >= 5 mn (26) 



2. The special first-order operator 



D:=(x-x')^- (27) 



has the Taylor basis (|T6| ) as a complete system of eigenfunctions: 



D\n >= n\n > for n = 0,1,2, (28) 

As a conclusion, the matrix elements of D read 

< m\D\n >= nd~ mn ; (29) 

hence, 

oo 

D = n\n X n\. (30) 



n=0 
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More generally, for an arbitrary polynomial P we have 

oo 

P(D) = J2 p (n)\n >< n\. (31) 

n=0 

This operator can be inverted provided the polynomial does not have roots at the 
integer positive points, i.e. if P(n) ^ for n — 0, 1, 2, . . . then 

OO 1 

W^E^X"! (32) 

The polynomial values P(n) in the denominator even improve the convergence if 
P(D) -1 is applied to an analytic function. 

3. The one- dimensional Schrddinger operator 

d 2 

L:=-j-£ + u, u — u(x), (33) 

has the matrix elements 

(Tfi \ 
\u m _ n for n < m, (34) 

< m\L\m + 1 >= 0, (35) 

< m\L\m + 2 >= 1, (36) 

< m\L\n >= for n > m + 3, (37) 

where 

"» : = (38 > 

The result can also be compactly expressed by: 

< m\L\n >= I j-u m _ n + < m + 2\n >, (39) 



if we set 

/ m 



:= 0, u n := for n < 0. (40) 
n I 
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It is also helpful to arrange the matrix elements to an infinite matrix 

/ 



' < m\L\n > 



w 
w a 
u 2 

U3 
U I 

a 5 





w 

Ml 



1 



W 



C) 




W 2 
ill 



G) 
G) 
G) 




1 


w 

G) 
G) 



U-2 






1 



w 



G) 



Wi 







1 



w 








1 





«« (;)«b (;)«8 (j)«a (s)«i «o 



v 



(41) 



/ 



This matrix can be expanded as an infinite sum of constant matrices with shifted 
diagonals multiplied by derivatives of u: 



1-2 + u I + uxh + ■ 



2 J i-2 



j>0 



where we formally set 

w_ 2 := 1, W-i := 0, 
and the proper meaning of each matrix Ij is just defined by ([4~1|). 

4. More generally, let us consider a linear differential operator of the form 



(42) 
(43) 



dx 1 



+ w, 



u = u(x). 



(44) 



where r is a positive integer. We will call it a higher-order Schrddinger-like operator. 
The one- dimensional Schrodinger operator is the special case r = 2. The calculation 
of the matrix elements in the preceding item can be extended to general r. We obtain 



< m\L r \n >- 



m 
v n , 



u, 



,+ <m + r n >, 



if we adopt the convention ( fH]|) . This defines the matrix representation of L r 

L r =< m\L r \n >= w 3 _ r J 3 _ r 

where now we alternatively set 

w_ r = 1, w_ r _i = ■ ■ • = w_i = 0, 
and the matrix J_ r is defined by 

< m\I- r \n >=< m + r\n > . 



(45) 
(46) 

(47) 
(48) 
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5. Finally, let us consider the powers L k of the Schrodinger-like operator (f44|) . Clearly, 
the matrix elements < m\L k \n > of the powers of the operator L r are given just by 
the k-th power of the infinite matrix L r (|46|). Thus 

L k = \m >< m\L k \n >< n\, (49) 

m,n>0 

where 

<m\L k \n>= Hi) u h-r ■ ■ ■ Uj k - r . (50) 

Here j = (ji, ■ ■ ■ ,jk) is a multiindex formed from non-negative integers j\, j 2 , ■ ■ ., 
jk and the constants /i(j) are determined by the products of the matrices If 

h(j) =< m\I h - r ■ ■ ■ Ij k - r \n > (51) 

3 Derivation of the KdV polynomials 

Let us apply a special multi- index formalism in order to present our result: we set 

m := (m 1 ,m 2 , . . .,m d ), (52) 

|m| p := mi + m 2 + • • • + m p , for p = 1,2, ... ,d, (53) 
\m\ \= \m\ d = mi + m 2 ^ h m d , (54) 

Theorem 1 The Taylor coefficients of the heat kernel coefficients = ak(x,x') to the 
one- dimensional Schrddinger operator (W read 

k 

< n\a k >= Y c ( 

(55) 

d=l \m\=n+2k-2d 

where the innner sum runs over integers ni\ > 0, m 2 > 0, . . . , m d > such that 

|m| = rrii + m 2 + • • • + fn d = n + 2k — 2d (56) 
and where the numerical coefficients read 
c(m) = c{m 1 ,m 2 , ■ ■ ■ ,m d ) 



The sum in ( \5lj) runs over integers i\, i 2 , . . . , i d _\ such that 

= i < ii < i 2 < ■ ■ ■ < id-i <i d = k, (58) 
Zip < |m| p+ i + 2p. (59) 
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Proof. The recursive differential equations ([H]) have the symbolic solution 



1 



1 



-i 



D L 



1 



1 



D 



-i 



L i l + jD 



-i 



L ■ 1 



(60) 



k ) V ' 2' 

This can be translated, by use of the eq. (|3~2"D, into a formula in terms of matrix elements: 



1 + i 



< n\a k >= 

x < nJLlrifc-i > • • ■ < n2\L\ni >< ni|L|0 > 



1 + T)"( 1 + f" 



(61) 



where nj, = n and the summation indices n\, n 2 , . . . , n&_i run through 0, 1,2,.. .. 

Recall the values (|34D -(|37D of the matrix elements of L. Although the sum in (|6ll) is 



formally an infinite one, it is effectively finite. Let us call a matrix element < m\L\n > 
essential if n < m and non-essential if n > m + 1. The number of essential elements in a 
term of (^) equals the degree d of the differential monomial in u. Let us analyze a typical 
term of degree d = N + 1 in (j3l|). Clearly, < ni|L|0 > is essential since n\ > 0. Let 



< n iN+1 \L\n iN >, . . . , < n il+ i\L\ni >, < n h \L\0 > . 



(62) 



denote all the essential elements, numerated from the right to the left. The space be- 
tween two essential elements is filled by a chain of unessential matrix elements, which can 
be written as Kronecker symbols. The unessential chain between < n^+^Lln^ > and 
< ni|L|0 > reads 

< n h + 2\n h > ■■■ <n 3 + 2\n 2 >< n 2 + 2\m > . (63) 

This product of Kronecker symbols cancels the summations over n 2 , n 3 , . . . ,n i± by fixing 
these running indices: 

n 2 = n x - 2, 

^3 = n>2 — 2 — m — 4, 



= n h _i - 2 = ... = n 1 - 2(*i - 1). 
The chain under consideration produces the numerical coefficient 



Co 



'i 
k 



1 1 



By virtue of fl64|) , it becomes 



1 + 



1 



1 + 



1 



n ix + i\ n 3 + 3 n 2 + 2 n x + 1 



Co 



f n x + r 



(64) 



(65) 



(66) 
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The next unessential chain between < ni 2+ i\L\ni 2 > and < ni 1+ i\L\ni 1 > reads 

< n i2 + 2|n i2 _i > • ■ ■ < n il+3 + 2\n il+2 >< n il+2 + 2\n il+1 > 
It cancels the summation over . . . , rii 2 by fixing 

nil+2 = ^ii+l — 2, 

n il+3 = n h+2 - 2 = n il+1 - 4, 

n i2 = n i2 _ x - 2 = ...= n h+1 - 2(i 2 - h - 1) 
and produces the numerical coefficient 



%2 / V ii+3/ V 2i + 2/ V zi + 1 
22 3 i\ + 2 + 1 



n i2 + z 2 n il+3 + ii + 3 n il+2 + Zi + 2 n n+1 + l 
By virtue of (|68|) , the latter becomes 



i 2 \ +*i + 1 



The procedure carries on this way; the last unessential chain fixes 

n i N +2 = n i N +l — 2, 

n iN+3 = n iN+2 -2 = n iN+1 - 4, 

n fc _i = n k -2 - 2 = . . . = n iN+1 - 2(k -i N -2) 
n k = n = . . . = n iN+ i - 2(k - i N - 1), 

and produces the numerical coefficient 

\in)\ k-i N J 

After the reduction procedure the typical term in (|61]) reads 

• • -cic < n iN+1 \L\n iN > ■■■ < n il+1 \L\ni >< ni|L|0 > 

and the summation indices are now N, i%, i 2 , iN, n i, n h+i, n n+\i ■ ■ ■ 

Notice, further, that n k = n is fixed. 
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In order to get a better expression for the essential matrix elements it is convenient to 
pass to new summation indices N, i±, i 2 , • • • , %n, m i, m 2, ■ ■ ■ , m N+i, where 

mi := Hi, 

m 2 := n il+1 - n h = n il+1 - m + 2{% x - 1), 

™3 := n h+1 - n i2 = n i2+ i - n ix+1 + 2(i 2 -i\ - 1), 

mAT+i := ra ijv+ i - n iN = n iN+1 - n iN _ 1+1 + 2(i N - i N -i - 1). (74) 

The old summation indices are expressed in terms of the new ones by 

ri\ = mi, 
n il+1 = mi + m 2 - 2i x + 2, 
n i2+l = m x + m 2 + m 3 - 2i 2 + 4, 



n iN+l = mi + m 2 -\ h m N+ i - 2i N + 2N. (75) 



Hence 



c 



Cl 



V + A 1 

'i 2 \ fmi + m 2 — i 3 + 3\ 1 



fk\fmi + --- + m N+1 -i N + 2N+l\ 1 

Some restriction for m-i, m 2 , . . . , m.jv + i is to be considered: 

mi + m 2 + --- + m N+ i = n + 2(k - N - 1); (77) 
it follows from (|74[) . Besides, all m^ are non- negative and we have to require 

< 2i p < mi + m 2 + • • • + m p+ i + 2p. (78) 
Finally, we express the essential matrix elements in terms of mi, m 2 , . . ., m^+i' 

( m i\ 

< m|L|o > = 1 iu mi , 

/mi + m 2 - 2«i + 2\ 
< n il+1 \L\n h > = \u m2 , 

\ m 2 



mi + m 2 + m 3 - 2z 2 + 4 

™3 



(79) 

in ^ fm 1 + --- + m N+1 -2i N + 2N\ 

< n iiV+ i|L|n ijv > = hW+i- ( 80 ) 
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We complete the proof by setting d := N + 1, «q := 0,id '■= k and summation over 
d — 1, 2, . . . , k. The asserted formula (|55| ) -(|57| ) follows. 

The outer summation in flSo] ) is a decomposition into homogeneous parts. Note, fur- 
ther, that the differential polynomial < n\a k > in u is isobaric, that means every term 
has the same weight 

|m| = mi + m 2 + • • • + = n + 2(/c — d). (81) 

Theorem 1 is a very special one-dimensional case of the general formula for the Taylor 
coefficients of the heat kernel coefficients to an arbitrary second-order partial differential 
operator of Laplace type obtained by one of the authors (I.G.A.) in || ||. We present 
here for the first time the complete proof; in |3J] the combinatorial details were omitted. 

Specializing to 

< 0\a k >= a k (x,x), (82) 

and considering (|12|) we obtain 

Theorem 2 The right-hand sides of the KdV hierarchy (Tjjj are given by 




where the multiindex m and the numerical coefficients c(m) have the same meaning as in 
Theorem 1. 

Notice that (|53p is valid for the matrix KdV hierarchy, since our proof did not make 
use of a commutative law. The convential scalar KdV hierarchy, where commutativity is 
assumed, emerges as a special case. 

4 On the structure of the KdV polynomials 

The explicit expression (^) is well suited to be evaluated by means of computer algebra. 
It is not practical to calculate higher terms by hand. Also, ([83]) conceals that some terms 
in Gk[u] have a relatively simple structure, as we are going to discuss now. The asymptotic 
expansion of the heat kernel diagonal 

1 

K(s; x, x) ~ (4^)-^ J2 tM x ' x ) sk = ( 4 ^r V2 2 £ 77^ G M s * ( 84 ) 

is a generating function for the sequence of the KdV polynomials G k - It is well known, 
that the diagonal value K(s; x, x) of the of the heat kernel obeys 



K(s;x,x) = e SU0 K(s;x,x), 



(85) 
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where K is the heat kernel diagonal calculated for u Q = 0. As a consequence, there holds 

\u l a k -i, (86) 



and 

where 
Notice that 



S9) 



G k = iz{ k (tooYGk-t, (87) 

i=0 V * / 

Ofe = afc| Uo= o, G k = G k\ Uo=0 - ( 88 ) 

V i J i\r(k-i + \) 2*\2i)\i)\i) 
Herefrom we find the terms of highest degree and lowest order in G k as far as we want: 

60 ^ G k = 30ug + 10 Q u k - 2 u 2 + 3 ^ u*- 3 (u 4 + 5u?) + • • • (90) 

The terms of lowest degree and highest order are found by means of some other general 
relation. Namely, we have shown in || that for the matrix KdV hierarchy the recursion 

k 

G fc = 2^Z 2i _ 1 G fc _ i , (k — 1,2,.. .), (91) 

i=i 

holds, where the differential polynomials Z k = Z k [u] on their side are recursively defined 
by 

Zk+i '■= -rZ k + 2. ZiZk-i, (92) 
dx £J 

Z x := u , (93) 

(The proof is far from being trivial in the non-commutative case.) This enabled us in |3[] 
to find 

G k = u 2k _ 2 + | y • + 1 J + (-!)*! «i«2fc-4-< + • • • , (94) 

where the points indicate terms of higher degree. 

Let us add here another result on the structure of the KdV hierarchy: we find, in the 
commutative case, all monomials in G k which are built solely from uq, U\ and u 2 . 

The proof uses again diagonal values of the heat kernel as generating function and the 
fact that Mo, U\ u 2 commute with each other. Besides, we will use the Bernoulli numbers 
B n , {n = 0, 1, 2, . . .), (cf., e.g., fl5|]) defined recursively by 

E^^^O, B :=l. (95) 

The first Bernoulli numbers read: B% = — |, B 2 = — |, B 3 = 0, -E> 4 = — ^. Generally, 
B 2k+1 = 0, B 2k = (-l)^ 1 !^! ^ for k > 1. 
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Theorem 3 For the scalar KdV hierarchy there holds 

where the inner sum runs over integer mi > 1, . . . , > 1 such that |m| = mi + • ■ ■ + 
rriN = k, the monomials bk are defined by 

h = u , (97) 

b 2n = -2 3n ~ 2 \B 2n \ul (98) 
n 

b 2n -i = -2 n -\2 2n ~l)\B 2n \u n 2 - 2 ul (n>2), (99) 
n 

and the dots indicate terms which effectively contain higher derivatives U3, etc. of u. 

Proof. Actually, we have to evaluate the asymptotic expansion of the heat kernel for a 
quadratic potential (harmonic oscillator) 

u(x) = u {x') + ux{x')(x - x) + ^u 2 (x'){x - x'f. (100) 

The heat kernel diagonal in this particular case reads (up to exponentially small terms 
that do not contribute to the generating function) 

K(s; x, x) ~ (4?rs)~ 1/2 exp jsu + $(s) + is 3 M^(s)| (101) 



where 



^ / x 1, /sinh(2su;)\ ,„„„x 

$ s = --log —± >- , (102) 

2 y Isuj j 

9(s) = SUJ -^ {SU) (103) 



and 



u = ^-u 2 /2 (104) 

Actually, this is the well-known heat kernel diagonal of a the one-dimensional harmonic 
oscillator, which can be read from the literature (cf., e.g. |TT| , Expanding the heat 
kernel diagonal one finds all the coefficients ak{x, x) in the considered approximation. The 
Taylor expansions of the functions $(£) and ^{t) read 

o3n— 2 1 f) I 

w-xhff**"- (105) 
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vis) = e ( , 9 ,; mn] ur 2 s 2n -\ (ice) 

n>2 \ An )- 

where B n are the Bernoulli numbers. Obviously, these series have only positive coefficients. 
We obtain the heat kernel diagonal in the form 

K(s; x, x) ~ (4ns)- 1 / 2 exp f E J , (107) 

where the coefficients b n are defined by (0)-(0). Note that all b n are positive. The 
heat kernel coefficients ak{x,x) and, therefore, the differential polynomials ( p4|) are 
determined by the expansion of the exponent, which completes the proof. 

Let us note that the heat kernel coefficients ak{x,x), in the above approximation, are 
given by the so-called Bell polynomials Bk(yi, ■ ■ ■ ,yk)- Namely, the Bell polynomials are 
defined by their generating function 

oo k / oo n \ 

ETF** = «P E^n ■ (108) 
k=0 K - \n=l m J 

Therefore, the formula (^BJ) can be also rewritten in the form 

1 (2k\ 

Gk[u] = -{ k jB k (b 1 ,...,b k ). (109) 
Several explicit expressions for the Bell polynomials are known, in particular, 

B n ( yi ,...,y d ) = ± E rif^" 1 )^ (U0) 

d=l |m|=nP=l V I" 1 !?- 1 / 

for n > 1 | I7| . Here, as usual, m = (mi, . . . , m^) with non- negative integers mi, . . . , 
and 

|m| p := mi + • • • + rn P , for p = 1, . . . , d (HI) 
|m| := 0, |m| := \m\ d (112) 



The expression ( |110| ) is valid for generally non-commuting variables yi,y2, ■ ■ •• Another 
expression, which holds only in the commutative case, is due to Faa di Bruno (cf. e.g. 
0): 

B n (yi,...,y n )=n\ E ft — f ' (H3) 

||m||=ni=l m i \3- / 

where the sum runs over all non- negative integers mi, . . . , m n such that 

| |m| | := mi + 2m-2 + . . . nm n = n (114) 

One should point out here that the theorem 3. is a particular case of a general result 
|16| of one of the authors (I.G.A.). In that paper the heat kernel diagonal to a Laplace 
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type operator L = — g MZ/ V^Vjf + Q acting on sections of a vector bundle V over a flat 



n- 



dimensional manifold with a metric g and the bundle connection V y was studied and the 
part of the heat kernel diagonal that depends solely on the curvature TZ^ = [V M , V„] of 
the bundle connection, the potential Q and its first V^Q and second V ' U Q derivatives 
is calculated in a closed form. Note that there is a misprint in the eq. (3.31) of the paper 

m 



5 Powers of Schrodinger-like operators 

The matrix elements of powers of a Schrodinger-like operator are already given in compact 
form by fl5"0|). Let us present here a more explicit expression for the powers of the proper 
Schrodinger operator (j33"l) . 



Theorem 4 The matrix elements of the k-th power of the r-th order Schrodinger-like 
operator ffi^l ) read 

k 

< m\L k r \n>=< m + rk\n> + E E K,k( m ) u m d • • • u mi (115) 

d=l m|=m— n+r(k— d) 

where the inner sum here runs over non-negative integers mi > 0, . . . , rrid > 0, such that 

|m| = m\ + • ■ ■ + rrid = m — n + r(k — d) (H6) 

and the numerical coefficients read 

7 / x_, / \ A (\m\ p + n-r(i p -p+l)\ 

h r>k (m) = h r) k{m 1 ,m 2 , . . . ,m d ) = 2^ 11 I • ( 117 ) 

i u -,i d ;p=l \ m P I 



The sum in runs over non-negative integers ii,...,id such that 

< h < ii < ■ ■ ■ < i d -i <i d <k, (118) 

|m| p + n - r(i p -p + 1) > 0. (119) 

The proof is performed by the method of matrix elements, analogous to the proof of 
the Theorem 1 and, hence, is omitted. 

Let us translate the result into the conventional notation. 



Corollary 1 The k-th power of the r-th order Schrodinger-like operator ( \^\ ) applied to 
a function f gives 

k 

L k r f = frk + E E h /mo (120) 

d=l |m|=r(fc-d) 
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where f m := rrzzf, m = (^o? m ) = ( m o, m i, • • • , ^d), the inner sum here runs over 
non-negative integers itlq > 0, ni\ > 0, . . . , > ; such that 

|m| = mo + mi + • • • + m^ = r(k — d) (121) 

and the numerical coefficients read 

h rik (m)=h rtk {m ,m 1 ,m 2 ,...,m d )= E II ( '"^ * %V P+1 j- ( 122 ) 

h,-,i d ;p=l \ m P ) 



The sum in ( \12$\ ) is the same as in the theorem 4- 



Proof. There holds, with a mixture of matrix element notation and conventional nota- 
tion, 

L k r f= E <0\L k .\n> f n . (123) 

0<n<rfc 

Here we insert < 0|L^|ra > from the theorem 4, i.e. we set m = there, rename mo := n 
and define a new multiindex m = (m , m). 

The numerical coefficients in ( |1 2U| ) for d — 1 and for d = k exhibit a simple structure, 
namely, 

k - x ( ri\ 

ftr,i(mo,mi) = E > ( i24 ) 

i=o \ m i/ 

h r ,k(m Q , mi, ... , m k ) = h r , k (0, . . . , 0) = 1. (125) 
Thus ^ 

L k J = frk+ E J2( m )u m J mo + --- + u k f (126) 

m + mi =r(fc-l) i=0 \ m l/ 

where the dots indicate the terms of degree in u greater than 1 and less that k. 

More explicit formulas for L k f are available for r — 1. The powers of the first order 
operator 

I*-±+n (127) 

are given by [[L7j 

L\f = jz{^B l {u ,... ) u l . l )f k . h (128) 

where B% = Bi(yi, . . . , yi) (I = 0, 1, 2, . . .) denotes the Bell polynomials, which we already 
discussed above. 
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